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Abstract 

The inviscid and thin accretion disc is a simple and well understood model system in accretion studies. 
In this work, modelling such a disc like a dynamical system, we analyse the nature of the fixed points of the 
stationary solutions of the flow. We show that of the two fixed points, one is a saddle and the other is a 
centre type point. We then demonstrate, using a simple but analogous mathematical model, that a temporal 
evolution of the flow is a very likely non-perturbative mechanism for the selection of an inflow solution that 
passes through the saddle type critical point. 
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1 Introduction 

In astrophysical fluid dynamics, studies in accretion processes occupy a position of prominence. Of such pro- 
cesses, critical (transonic) flows — flows which are regular through a critical point — are of great importance 
A classic example of such a flow, in the stationary spherically symmetric situation, is the Bondi flow '7. What 
is striking about the Bondi solution is that while it is virtually impossible to numerically generate it in the sta- 
tionary limit of the hydrodynamic flow, it is easy to lock on to the Bondi solution when the temporal evolution 
of the flow is followed. This dynamic and non-perturbative selection mechanism of the transonic solution is in 
very favourable agreement with Bondi's conjecture that it is the criterion of minimum energy that will make 
the transonic solution the favoured one |21ini- 

The appeal of the spherically symmetric flow, however, is somewhat pedagogic, limited as it is by its not 
accounting for the fact that in a reasonably realistic situation, the infalling matter would be in possession of 
angular momentum — and hence the process of infall should lead to the formation of what is known as an 
accretion disc. A well understood and quite regularly invoked model of such a system is the inviscid and thin 
accretion disc In this work, we show how such a physical system — studied in its stationary limit — 

gives rise to saddle point like behaviour for at least one of the critical points of the flow. The adverse implications 
regarding the realizability of a solution passing through such a point is then addressed by considering the system 
in a real non-linear time evolution. A simplified mathematical analog of such a situation shows that the evolution 
indeed selects the critical (transonic) solution. Finally, we discuss that in the vicinity of the critical point, the 
resulting perturbation equation arising out of a linear stability analysis of the stationary solutions, quite closely 
resembles the metric of an acoustic black hole, and that is probably the signal that trajectories pass through 
the critical point in a unidirectional fashion. 

2 The stationary equations of the flow and its critical points 

The flow system to be considered is an axisymmetric thin disc of a compressible astrophysical fluid, with a 
thickness _ff 7 . In the vertical direction the flow is considered to be in hydrostatic equilibrium. The governing 
equations for the flow should therefore be the equation of continuity and the equations for momentum balance 
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in both the radial and the azimuthal directions. It is customary to write the stationary flow equations in the 
following notation 0. 



Equation of continuity: 



• Radial momentum balance equation: 

dv 

• Angular momentum balance equation: 
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In the above, v is the radial velocity, R is the radial distance, il is the angular velocity, fix is the Keplarian 
angular velocity defined as flx^ = GM/R^, Cg is the velocity of sound defined as c^^ = dP/dp and a is the 
effective viscosity of Shakura and Sunyaev [S] . 

It is a standard practice to make use of a general polytropic equation of state P = kp^ where k and 7 are 
constants, with 7 being the polytropic exponent ^U], whose admissible range (1 < 7 < 5/3) is restricted by the 
isothermal limit and the adiabatic limit respectively. The condition of hydrostatic equilibrium along a direction 
perpendicular to the thin disc, allows for the use of the approximation {H / R)={cs/vk), where vk = R^k El- 
In that case Eq. (1) leads to 



"^wi?^/^ = constant 



where n = (7 — 1) ^ . Combining Eqs. (2) and (4) yields, 
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and a first integral of Eq. (3) can be written as 
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where L is a constant of integration. 

We now consider the inviscid limit of Eq. (6) by setting a = 0. This, somewhat simplifying, prescription has 
found regular favour in accretion literature 0)13112], and it gives the condition UR^ = L, where L, which can 
be physically identified as the specific angular momentum of the flow, now becomes a constant of the motion. 
Such a constraint allows us to fix the critical points of the flow. At such points both the numerator and the 
denominator in the right hand side of Eq. (5) vanish simultaneously ^E], and this will deliver the critical point 
conditions as 
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with the subscripted label c indicating critical (fixed) point values. We should be able to find Vc and Rc from 
Eqs. (7). The second relation above, which is actually a quadratic in Rc, leads us to 
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and this gives the two critical points at radii R — Rci and R = Rc2, with i?c2 > Rci- To fix the critical point 
coordinates {Rc,Vc), we need to look at the integral of Eq. (2), which reads as 
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Figure 1: A schematic diagram of the accretion problem. The fixed point C2 is a saddle point. Ci is a centre type 
point. The separatrices pass through C2. 



in which the constant of integration has been determined with the help of the boundary condition that for very 
great radial distances, the speed of sound approaches a constant "ambient" value Cs (00), while v — > 0. In this 
expression wc use the condition Vc'^ = 2csc^/(7 + 1) the value of Rc from Eq. (8), to fix Csc in terms of the 
constants of the system. This will in turn fix the critical point coordinates {Rc,Vc)- 



3 The disc as a dynamical system 



To have any understanding of the nature of the critical points of the inviscid and thin disc being studied here, 
we need to cast the equation giving the stationary solutions of the disc system, along the lines of a dynamical 
system. To do that we appeal to Eq. (5) and parametrize it in the form 
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It is to be noted that the parametrization has been carried out in an arbitrary mathematical parameter space 

r, and in such a space the two equations above represent a dynamical system. To analyse the nature of the 
fixed points, we expand and linearize about the fixed point coordinates (vc^jRc) in Eqs. (10), and in terms of 
the perturbed quantities we obtain a set of linear equations given by 
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Use of solutions of the form e^'^ would deliver the eigenvalues A — growth rates of Sv^ and SR — as 
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where /3 = (5 — 37) /(5 — 7) and the local Keplerian angular momentum is Lkc = ^/GMRc- It is now quite 
evident that if Rc is the outer fixed point Rc2, then A^ is positive and the fixed point is saddle type, while if Rc 
is the inner point Rd, then A^ is negative and the fixed point is centre type. We show the trajectories in Fig.l, 
in which the arrows indicate the direction of flows. 

The fixed point C2 at i? = Rc2 is the saddle. In a manner of speaking, it may be dubbed the "sonic point" 
because one of the two trajectories passing through it is transonic, rising from subsonic values far from it to 
supersonic values when R < Rc2- This trajectory is shown as the heavy solid curve, and is the accretion flow. 
The heavy dashed curve is the wind. The above demonstration that the fixed point C2 is a saddle has adverse 
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implications regarding the realizability of the transonic trajectory passing through it. Conventional wisdom 
about the nature of saddle type points [HI gives us to believe that if we go directly to the static limit of the 
dynamic evolution, then the transonic trajectory would not be realizable. However, it is generally believed that 
transonic trajectories do occur in astrophysical systems such as the one being studied here jlj. Our contention is 
that it is the temporal evolution of such a system that selects the transonic trajectory. This selection mechanism 
is purely non-perturbative as opposed to the perturbative technique of a linear stability analysis of all physically 
feasible stationary inflow solutions. The latter method offers us no clue as to the choice of the system for any 
particular solution, since all stationary solutions have been found to be stable under the influence of a linearized 
perturbation — considered both as a travelling wave and a standing wave jl^' . 

As an aside it may also be mentioned here that the nature of the outer fixed point C2 in Fig.l, has a bearing 
on the possible range of values that the constant specific angular momentum L may be allowed. For solutions 
passing through the outer critical point, it can be well recognized that the flow would be sub-Keplerian 00]. 
In that case, it may be concluded that the condition < {L/Lkc2Y < 1 would hold good. In addition to that, 
from Eq. (12), the saddle type behaviour of the point C2 would also imply that there would be another upper 
bound on L, given by (LjLKc'i)^ < /3. For the admissible range of the polytropic index 7, the possible range 
for /3 would be < /3 < 1/2. This would then imply that the latter bound on L would be more restrictive, as 
compared to the former. It is interesting that this essentially physical conclusion could be drawn from modelling 
the disc as a dynamical system in a mathematical parameter space. 



4 A model for a dynamic selection mechanism 

In the previous section we contended that the temporal evolution of the disc system selects a critical solution 
that passes through the outer critical point. To understand how this comes about, we consider a simplified 
mathematical example, since it is well known that the equations of a compressible flow cannot be exactly 
integrated. 

The model system that we introduce, describes the dynamics of y(x^ t) as 

^ + iy-x)^=y + 2x + x^ (13) 

whose static limit yields 



dy y + 2x + x'^ 
dx y — X 

and which, viewed as dynamical system, is seen as 



(14) 
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y + 2x + x'^ 

y-x (15) 



In the y — X space, the flxed points {xc,yc) are to be found at (0,0) and (—3,-3). As in the preceding 
section, a linear stability analysis of the fixed points in r space gives the eigenvalues A by = 1 + 2{xc + 1)- 
It is then easy to see that (0, 0) is a saddle point while (—3, —3) is a centre type point. The integral curves are 

2 

- 2a;?; - 2x'^ - -x^ = c (16) 

and the trajectories passing through the saddle point are the ones with c = 0. The different possible trajectories 
are shown in Fig. 2. The separatrices are shown as the heavy solid curve and the heavy dashed curve. The 
similarity between Figs.l & 2 is obvious. We now explore the dynamics. 

To obtain a solution to Eq. (13), we need to apply the method of characteristics ^21- This involves writing 

(U ^ dx ^ dy 
1 y — X y + 2x + x'^ 

The task is to find two constants ci and C2 from the above set and the general solution of Eq. (13) would be 
given by ci = _F(c2), where the function F is to be determined from the initial conditions. It is easy to see that 
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Figure 2: A schematic diagram of the model problem. The curve given by Eq. (16) passes through the origin for c = 0. 
For this critical trajectory O is a saddle point while C is a centre type fixed point. This figure differs from Fig.l by a 
tilt of the axes. 



one of the constants of integration is clearly the c of Eq. (16). Hence, we write ci = c and from the first part 
of Eq. (17), we can write 

dt = ±[ '^'^ = (18) 

J ^3x2 + (2/3)a;3 +c 

which solves the problem in principle. To put this in a usable form, we need to carry out the integration of 
Eq. (18). This cannot be done exactly. For small x (the most important region, since it is near the saddle), 
we can drop the x^ term to a good approximation. Further, we choose only the positive sign in the right hand 
side of Eq. (18) by the physical argument that we would wish to evolve the system through a positive range of 
t (time) values. Integration of Eq. (18) will then lead us to the result 
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which will then make the solution of Eq. (13) look like 
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We wish to study the evolution of the system from the initial condition y = for all x at t — 0. For small x, 
dropping the x^ term again, allows us to determine the form of the function F as F{z) ~ —3(2 
solution, consequently, becomes 
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and as t — >cxd, we select the separatrices given by 




(21) 
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(22) 



a result whose remarkable feature is worth stressing. The evolution started far away from the separatrix. In 
fact, it started as y = for all x &t t — 0. The evolution proceeded through a myriad of stationary states 
(all arguably stable under a linear stability analysis) and then locked on to the separatrix. Our claim is that 
this dynamic mechanism selects the velocity profile in the practical situation as well, and this is entirely in 
consonance with a similar and well established selection mechanism for the spherically symmetric flow [H]. 
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5 Concluding remarks 



In a previous section we discussed that linear stability analysis of the stationary inflow solutions offers no direct 
clue as to the preference of the disc system for a particular stationary solution. However, in a somewhat indirect 
approach to this issue we could make use of a linear stability analysis to have a feel for the special feature of the 
critical solution. A very convenient way of carrying out the linear stability analysis is in terms of the variable 
<p = p^^+^^/^wi?^/^ [21 • In terms of this variable, the perturbation equation for reads 



a<2 dR\ ' dt va dR 



= (23) 



in which (j)' is the perturbation on the steady solution of 0, and the subscripted label indicates steady state 
values wherever in use [T^ . 

In a spacetime with metric g^^'^ , the wave equation for a scalar variable is ^1] 

Ignoring the angle variables, and using the notation that /i = 1 corresponds to t and /i = 2 corresponds to R, 
we have the identification 



7 + 1 



The effective covariant components are 



1 vo 1 

311 = 1 - -r^, 9i2 = 921 = -r^, 922 = — — (26) 

CsO CsO CsO 



in which c~so = \/2(7 + 1) ^/^Cso- 

The differential distance in the metric space has the form 

ds^ = [cSo^df - [dx' - vo'dt) <5y {dx^ - vo^dt)] (27) 

CsO 

and this is to be compared with the Painleve-GuUstrand form of ds"^ for a rotating black hole of mass M 

ds^ = _ I ^^^^1^^ -^dt] + r^dO^ (28) 

We now note from Eq. (7) that at the sonic point, Vc^ = {2/b)[{GM / Rc) — {Lj R^Y], and hence in the vicinity 
of the sonic point the metric implied by the linear stability analysis of the stationary solution appears quite 
close in form to the metric of a rotating black hole. The unidirectional nature of trajectories near a black hole 
suggests that trajectories go through the sonic point, which is equivalent to saying that transonic trajectories 
will be chosen. 

In concluding we would also like to comment on the nature of the possible critical point (s) when we consider 
the more likely case of a disc with a mechanism for the outward transport of angular momentum. It is possible 
with the help of Eq. (6) to expand and linearize about the critical point value of the angular velocity 17, but 
without any knowledge of an exact dependence of f2 on i?, there would be no way of fixing the critical point (s). 
However it can be conceived that as an extension of the inviscid case, there may be multiple parameter-dependent 
critical points and in that event at least one of them may be saddle type, with all its associated difficulties 
within the stationary framework itself. In that event, once again the selection mechanism has to be time 
evolutionary in nature. 
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